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Abstract 

We present an analytic example of two dimensional quantum mechanical system, 
where the exponential suppression of the probability of over-barrier reflection changes 
non-monotonically with energy. The suppression is minimal at certain "optimal" ener- 
gies where reflection occurs with exponentially larger probability than at other energies. 

1 Introduction 

Tunneling and over-barrier reflection are the characteristic non-perturbative phenomena in 
quantum mechanics. They typically occur with exponentially small probabilities, 

V oc e~ F/h , (1) 

where F is the suppression exponent; still, the above phenomena are indispensable in under- 
standing a wide variety of physical situations, from the generation of baryon number asym- 
metry in the early Universe [1] to chemical reactions [2] and atom ionization processes [3]. 

During the last decades extensive investigations of tunneling processes in systems with 
many degrees of freedom have been performed [2J, SI El El El El El OH EH EGJ [13] . These studies 
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revealed a rich variety of features of multidimensional tunneling which are in striking contrast 
to the properties of one-dimensional tunneling and over-barrier reflection. In particular, the 
following phenomenon has been observed: the probability of tunneling may depend non- 
monotonically on the total energy of the system and exhibit resonance-like peaks. One can 
envisage three physically different mechanisms of this phenomenon. The first mechanism, 
present already in one-dimensional case, is tunneling via creation of a metastable state. 
In this case the tunneling probability at the maximum of the resonance is exponentially 
higher than at other energies. On the other hand, the resonance width AE is exponentially 
suppressed; so, after averaging with an energy distribution of a finite width the effect of the 
resonance is washed out in the semiclassical limit H — > 0. The second possible mechanism of 
non-monotonic behavior oiV(E) is quantum interference [3 [13] (see also |14]). In this case 
the peak value of the tunneling probability is only by a factor of order one higher than the 
average value, while the width of the resonance scales as AE cx h. Again, the resonances 
become indiscernible in the semiclassical limit. In both these cases the resonances can 
be attributed to the subleading semiclassical corrections, i.e. non-monotonic behavior of 
the pre-exponential factor omitted in Eq. ([T]). The third possibility is that the suppression 
exponent F(E) is non-monotonic. In this case the existence of the "resonances" is the leading 
semiclassical effect: the optimal tunneling probability at the maximum of the resonance is 
exponentially higher than the probability at other energies. At the same time the resonance 
width scales asj AE oc y/h. This last possibility of "optimal tunneling" is definitely of 
interest; yet, it did not receive much attention in literature. We are aware of only a few works 
mentioning non-monotonic dependence of the suppression exponent on energy [151 ESI E]. 
It is worthwhile studying this phenomenon in detail; this can provide a new insight into the 
dynamics of multidimensional tunneling. 

In this paper we consider the process of over-barrier reflection in a simple model with two 
degrees of freedom. Our setup is interesting in two respects. First, the model under study 
is essentially non-linear and the variables cannot be separated; still, over-barrier reflections 
in this model can be described analytically within the semiclassical framework. Thus, this 
model can serve as an analytic laboratory for the study of multidimensional tunneling. Sec- 
ond, the suppression exponent F of the reflection process behaves non-monotonically as the 
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total energy E changes. We demonstrate that the function F(E) possesses a number of 
local minima E = E Q , where reflection is optimal. We stress that the process we study is 
exponentially preferable at "optimal" energies as compared to other energies. 

Our model describes the motion of a quantum particle in the two dimensional harmonic 
waveguide (see Refs. [BHHHEIl] for similar models). The Hamiltonian is 



where x, y are the Cartesian coordinates and m is the mass of the particle. The function 
U = muj 2 w 2 /2 represents the waveguide potential in two dimensions: a particle with small 
energy is bound to move along the line w(x, y) ~ 0. We do not introduce a potential barrier 
across the waveguide and consider the case when the line w = stretches all the way from 
x — > — oo to x — > +00. We also assume that the function w(x, y) is linear in the initial 
asymptotic region, 



In the present paper we consider two particular cases of the function w(x,y) describing 
waveguides with one and two sharp turng§, see Fig. [TJ 

The motion of the particle at x — > —00 is a superposition of free translatory motion in 
x direction and oscillations of frequency lu along y coordinate; the state of such a particle is 
fully characterized by two quantum numbers, the total energy E and ^/-oscillator excitation 
number N. The particle sent into the waveguide from the asymptotic region x — > — 00 with 
given E, N may either continue to move towards x — > +00, or reflect back into the region 
x — > —00. We are interested in the probability V(E, N) of reflection. 

Let us discuss reflections at the classical level. [Note that the classical counterpart of N 
is the energy of transverse oscillations.] Consider first the waveguide with one sharp turn 
(Fig. UK). One observes that the outcome of the classical evolution, i.e. whether or not 
the particle reflects from the turn, depends not only on the total energy E, but on other 
dynamical quantities as well. In particular, the direction of the momentum of the particle 
in the vicinity of the turn (point C on the graph) is important. This means that the entire 
dynamics in the waveguide should be taken into account in order to determine the possibility 
of classical reflection. This is in sharp contrast with the situation in one-dimensional case, 
where reflection from the potential barrier (or transition through it) is ensured by the value 
of the conserved energy of the particle. 

5 The explicit expressions for the waveguide functions w(x, y) will be presented in the subsequent sections. 
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Figure 1: The equipotential contour U = E for the waveguides with (a) one and (b) two 
sharp turns. An example of classical trajectory is shown in the case (b). 

Now, consider the waveguide with two turns. The model is characterized by the angles 
of the turns and the distance L between them (see Fig. [lb). Suppose the particle starts 
moving classically from x — > — oo with N = along the valley w — 0. Then, the transverse 
oscillations get excited only after the particle crosses the first turn, point C on the plot, so 
that at the time of arrival to the second turn (point C) approximately uot/2h oscillations 
are made, where r ~ LyJrn/2E is the time of motion between the two turns. The state of 
the particle (coordinates and momenta) at which it comes across the second turn depends 
periodically on the phase of transverse oscillations lot. Hence, one expects that the regime 
of motion of the classical particle can change from transmission to reflection and back as the 
energy grows (r decreases); the energies where it happens can be roughly estimated as 

muj 2 L 2 

En ~ 2(2^ ' (2) 

We will see that this is indeed the case for the waveguides with certain angles of the turns. 

At some values of E, N the reflection process cannot proceed classically. Then, at the 
quantum mechanical level its probability is exponentially suppressed, F(E, N) > 0. It is 
natural to call such a process "over-barrier reflection']^]. The central quantity to be studied 
below is the suppression exponent F(E, N) of this process. The above discussion suggests 
that F(E,N), being determined by the entire dynamics in the waveguide, may be a highly 
non-trivial function. For the particular case of the waveguide with alternating regimes of 



6 By this term we want to emphasize that the process is classically forbidden. Recall, however, that there 
is no actual potential barrier across the waveguide in our setup. 



4 



classical reflections and transmissions F should oscillate: F = at the energies where the 
classical reflections are allowed, and F > at the energies where the reflections are classically 
forbidden. One can expect that the similar oscillatory behavior of the suppression exponent 
persists for other two-turn models as well. Now, instead of reaching zero, F may possess 
a number of local positive minima implying that the reflection at the "optimal" energies is 
still a tunneling process. 

Let us emphasize the difference of the "optimal tunneling" from quantum interference 
and resonance phenomena in our two-turn model. The interference of the de Broglie waves 
reflected from the two turns can, in principle, lead to oscillations in the reflection proba- 
bility V(E). One can estimate the positions of the interference peaks by equating the De 
Broglie wavelength of the particle to an integer fraction of the distance between the turns, 
2nh/v2mE ~ L/n. This yields the energies of the interference peaks, 



This formula is completely different from Eq. ([2]) for the peaks due to "optimal tunneling" . 
In particular, the distance between the adjacent inteference peaks, 



scales proportional to h. Thus, these peaks should be averaged over in the semiclassical 
limit. Besides, the amplitude of the interference peaks is at most of order one and does 
not affect the suppression exponent. Indeed, the exponential increase of the scattering 
amplitude can arise due to quantum interference only in the presence of a resonant state 
with exponentially long life-time. This state should be supported somewhere in between the 
turns and should be classically stable. In Sec. 14.21 we show that such states are absent in 
our system. One concludes that the peak-like structure of the probability V(E) of "optimal 
tunneling" is caused by completely different physical reasons as compared to the case of 
resonance scattering in quantum theory. 

It is worth noting that the phenomenon of "optimal tunneling" has an important imple- 
mentation in field theory. Recently it was argued [T7j (see also Ref. [16]) that the probability 
of tunneling induced by particle collisions [T8j [19] reaches its maximum at a certain "op- 
timal" energy and stays constat at higher energies. This result, if generic, provides the 

As opposed to the quantum mechanical case, the tunneling probability does not decrease at energies 
higher than the "optimal" one. This is due to the possibility, specific to the field theoretical setup, to emit 
the excess of energy into a few hard particles, so that tunneling effectively occurs at the "optimal" energy. 
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answer to the long-standing question [20] about the high-energy behavior of the probability 
of collision-induced nonperturbative transitions in field theory. The quantum mechanical 
model presented here supports the generic nature of the phenomenon of "optimal tunnel- 
ing"; the simplicity of our model enables one to get an intuitive insight into the nature of 
this phenomenon. 

The paper is organized as follows. In Sec.[2]we review the semiclassical method of complex 
trajectories, which is exploited in the rest of the paper. Reflections in the waveguides with 
one and two turns are considered in Sees. [3] and H] respectively. We discuss our results in 
Sec. In appendix we analyze the validity of some assumptions made in the main body of 
the paper. 



2 The semiclassical method 



We start by describing the semiclassical methoqj of complex trajectories which will be used 
in the study of over-barrier reflections. We concentrate on the derivation of the formula 
for the suppression exponent F(E, N) (see Refs. [2J [8], [9] for the details of the method and 
Ref. [19] for the field theory formulation). In what follows we use the system of units 

h = m = uj = 1 , 

where the Hamiltonian takes the form, 

H=^(p 2 x +p 2 y + w 2 (x,y)) . (3) 

One starts with the amplitude of reflection into the state with definite coordinates 
x f <0, y f , 

A=(xf,y f \e- t *to-V\E,N). (4) 

Here \E, N) is the initial state of the particle moving in the asymptotic region X{ — > — oo 
with fixed translatory momentum p = \/2(E — N) and the oscillator excitation number N. 
Semiclassically, 

( Xi , yi \E, N) = ^ cos (J^_P y (y')dy' + rr/^j , (5) 



8 Note that the method has been confirmed by the explicit comparison with the exact quantum mechanical 
results in Refs. [H [14] ; specifically, the recent check [14] deals with the case when the dependence of the 
suppression exponent on energy is not monotonic. 
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where Xi, Hi denote initial coordinates, 

Py(y') = V2N-y", (6) 

and we omitted the pre-exponential factor which is irrelevant for our purposes. Using Eq. (jHJ), 
one rewrites the amplitude P| path integral, 

Xf,yf 

A = J dxdVi j[dx][dy] e iS+ipoXi cos (J y ^p y (y'W + it/ A , (7) 

where S is the classical action of the model fl3J). 

In the semiclassical case the integral (j7j) is dominated by the (generically complex) saddle 
point. Note that, as we continue the integrand in Eq. (j7j) into the plane of complex coor- 
dinates, one of the exponents constituting the initial oscillator wave function grows, while 
the other becomes negligibly small. Within the validity of our approximation, we omit the 
decaying exponent by writing 

cos ( / p y (y')dy' + ?r/4 ) -> exp (i ! p y (y')dy' 

\JV2N J \ J \f2N 

with the standard choice^] of the branch of the square root in Eq. ([6]). 

One proceeds by finding the saddle point for the integral (|7j) with the substitution (JHJ). 
Extremization with respect to x(t), y(t) leads to the classical equations of motion, 

'x = —ww x , y = —wWy . (9) 

Differentiating with respect to Xi = x{U), yi = y(U), one obtains, 

±i = po = ^/2(E - N) , yi =p y ( yi ) = ^2N -yl . 

The latter equations are equivalent to fixing the total energy E and initial oscillator energy 
iV of the complex trajectory, 

E = ^ + N, (10a) 
N = Uyf + yf) • (10b) 



9 The correct branch is fixed by drawing a cut between the oscillator turning points y — ±V2iV, and 
choosing lmp y > at y e E, y > V2N, see, e.g., Refs. [21] . 
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Substituting the saddle-point configuration^] into Eq. (J7|), one obtains the amplitude of the 
process with exponential accuracy, 

A OC e iS + iB ( x T.>Vi) 



where the term 



B{xi, jji) = p Xi + p y {y')dy' (11) 

JV2N 

is the initial-state contribution. For the inclusive reflection probability one writes, 

The integral over the final states can also be evaluated by the saddle point technique; ex- 
tremization with respect to Xf = x(tf), yf = y(tf) fixes the boundary conditions in the 
asymptotic future, 

Imxf = Imxf = , Imy/ = Imyj = . (12) 

In this way one obtains the expression ([1]) for the reflection probability, where the suppression 
exponent F is given by the value of the functional 

F(E, N) = 2 Im S + 2 Im B(x h y t ) 

evaluated on the saddle-point configuration — a complex trajectory satisfying the boundary 
value problem (j9|). (fTOjh (TT2l). 

The contribution B(xi, yt) of the initial state is simplified after one uses the asymptotic 
form of the solution at t — > — oo (xi — > — oo), 

x = Pot + x , y = ae~ lt + ae lt . (13) 

Equations (fTUl) guarantee that the quantities po = \J2(E — N) and 2aa = N are real, since 
E, N G M. Therefore, one may introduce two real parameters T, 9 as follows, 

2 Im x = -p T , a = a*e T+e . (14) 

One finds for the initial term f[TT|) . 



21m B(x h Vi ) = Im |2p z i - 2A^arccos(?/ i /v / 2Ar) + y iX /2N - yf 
-p 2 Q T-N(T + 9)+lm(y iyi ) , 



10 For simplicity we assume that the saddle-point configuration is unique. Otherwise, one should take the 
saddle point corresponding to the weakest exponential suppression. 
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and thus 

F = 21mS-ET-N6 , (15) 
where S is the classical action of the system ([3]) integrated by parts, 

1 f tf 

S=--j dt [xx + yy + w 2 (x,y)] . (16) 

Let us comment on the physical meaning of the parameters T,9. Consider two trajectories 
which are solutions to the boundary value problem (Q, (ITU1) . ffT2|) at neighbouring values of 
E, N. The differential of the quantity 2 Im S as one deforms one trajectory into the other is 

d (2 Im S) = d lm(2S + Xiii + yiy\) = lm(xidxi — x-idxi + yidiji — ijidyi) = EdT + Nd9 , 

where in the last equality we used the asymptotic form ffl3|) . (|14p of the solution. Then, 
from Eq. ([15"]) one finds, 

dF(E, N) = —TdE - 6dN . (17) 

Thus, the parameters T and 8 are (up to sign) the derivatives of the suppression exponent 
with respect to energy E and initial oscillator excitation number iV respectively. 

Our final remark is that the boundary value problem (Q, fflOl) . ffT2|) is invariant with 
respect to the trivial time translation symmetry, 

t^t + 5t, 5teR, (18) 

which can be fixed in any convenient way. 

3 The model with one turn 

To warm up, we consider the simplest model, where the waveguide has one sharp turn, 

w = y 9(—x + y tg (3) + cos f3 (x sin f3 + y cos f3) 9(x — y tg (3) . (19) 

Here 9(x) is the step function. It is convenient to use the rotated coordinate system, 

£ \ / cos j3 — sin (3 \ ( x 
7] J \ sin (3 cos (3 J \ y 

The waveguide function takes the form, 

w = r? cos sin /?#(-£) . (20) 
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The equipotential contour u> 2 (£, rj) = const is shown in Fig.[2j One observes that the motion 
of the particle in two regions, £ < and £ > 0, decomposes into the translatory motion and 
oscillations in the coordinates x, y and £, 77 respectively (see. Eqs. ( fl9]) and (120]) ): the 
frequency of //-oscillations in the latter case is cos (3. 

Due to the presence of the step function, the first derivatives of the potential fT2T)]) are 
discontinuous^"*] at £ = 0. Strictly speaking, the semiclassical method is not applicable in 



this situation [21]. Thus, the formula (120]) should be regarded as an approximation to some 
waveguide function with smooth turn. Generically the width of the smoothened turn is 
characterized by a parameter b; the sharp-turn approximation (1201) corresponds to b — > 0. 
An example of smoothening is provided by the following substitution in Eq. (1201) . 

m - m = Y+~pub ■ (21) 

The semiclassical description can be used as long as the de Broglie wavelength of the particle 
is small compared to the linear size of the potentia ^1,1/ \[E b. We conclude that the 
sharp-turn and semiclassical approximations are valid simultaneously for smooth waveguides 
with 

1>6>1/v A E. (22) 



n Note that the potential itself is continuous. 

12 Another semiclassical condition is that the energy is sufficient to excite a lot of oscillator levels, £> 1. 
It is satisfied provided Eq. ([22]) holds. 
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An important property of the model (|20|) is invariance of the classical equations of motion 
(IH]) under the rescaling of the coordinates, 

x -> Ax , y -> Ay . (23) 

Using the transformation (|2"3"1) . one may express a solution x(t), y(t) with energy in terms 
of the "normalized" one, 

x = xVE , y = y\fE , 
where the solution x(t), y(t) has unit energy; its initial oscillator excitation number is 

v = N/E . 

The suppression exponent ffl5l) takes the form, 

F(E, N) = Ef p (v) , (24) 

where fp{y) is the exponent for the "normalized" solution. Substituting the expression (jMj) 
into Eq. (ITT)) , one obtains, 

fp(v) = -T-6u. (25) 

We will exploit Eq. (|25|) in the end of this section. Now, we proceed to finding the "normal- 
ized" trajectories. 

At certain initial data v > u CT the particle can reflect from the turn classically, so that 

fp(u > p CT ) = . 

Let us find the value of z/ cr . In the region £ < the classical solution takes the form, 

x (t) = Pot + x , (26a) 
y{t) = A sin(t + <p) . (26b) 

Having crossed the line £ = (line AB in Fig. [2]), the classical particle can never return 
back into the region £ < 0. Indeed, in this case it moves at £ > with constant momentum 
> 0. Thus, the particle can reflect classically only if its trajectory touches the line 
£ = 0. The potential of our model has ill-defined derivatives at £ = 0, and the fate of the 
particle moving along the line AB depends on the particular choice of the smoothening of 
the potential. In appendix we consider the motion of the classical particle in the case when 
nonzero smoothening of width b is switched on. For a class of smoothenings we show that 
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in the small vicinity (5£ ~ b) of any trajectory touching the line £ = there exists some 
"smoothened" trajectory, which reflects classically from the turn. Consequently, below we 
associate the trajectories touching the line £ = with the classical reflected solutions. 
One notices that the inclination of the trajectory fl2B|) is bounded from above 



dy 



dx 



Po 



therefore, the classical trajectory of the particle can touch the line £ = 0, that is, y/x = ctg j3 
only at 

A)/S>ctg/3. (27) 

From Eqs. (|27|) . (|26[) . (fTUj) one extracts the condition for the particle to reflect classically 
from the turn, 

v > v cx = cos 2 (3 . (28) 

The critical classical solution at v — u CT touches the line £ = at rj = (point C in Fig. [2]), 
where its trajectory 

x a (t) = V2tsin(3 , (29) 
y C r(t) = v^2sintcos/3 . 

has the largest inclination. 

We now turn to the classically forbidden reflections at v < v CT , which are described by 
the boundary value problem (jHJ), ( |T0|) . (112p . One makes the following important observation. 
The waveguide function (120]) has the form of two analytic functions glued together at £ = 0. 
Hence, the equations of motion (Q can be continued analytically to the complex values of 
coordinates in two different ways, starting from the regions £ < and £ > respectively. In 
this way one obtains two complex solutions, £-(£), V-(t) an d C+(0> V+(t)- These solutions 
and their first derivatives should be matched at some moment of time t\, = 0. [Note 
that the matching time t% does not need to be real.] Below we conventionally refer to these 
solutions as the ones belonging to the regions £ < and £ > 0. 

By the same reasoning as above we find that once the particle arrives into the region 
£ > 0, it never reflects back to £ < 0, unless = 0. So, in the region £ > one writes, 

£+(t) = , (30a) 

r) + (t) = sin(tcos/3 + (p„) , (30b) 

cos p 
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where the "normalization" condition E = 1 has been used explicitly. Due to the conditions 
in the asymptotic future, Eqs. (TT2"j) . the parameter <p v is real. We use the translational 
invariance ffT8l) to set <p v = 0. Note that we again associate the trajectory going along the 
line £ = with the reflected one. 

The physical picture of over-barrier reflection that comes to mind matches with the new 
mechanism of multidimensional tunneling proposed recently in Refs. 0, [TT] . The process 
proceeds in two steps. The first step, which is exponentially suppressed, is formation of the 
periodic classical orbit (!30l) oscillating along the line £ = 0. This orbit is unstable. At the 
second step of the process the unstable orbit decays classically forming a trajectory going 
back to x — > — oo at t — » +oo. Clearly, the second step does not affect the suppression 
exponent of the whole process, and we do not consider it explicitly. In what follows we 
concentrate on the determination of the tunneling trajectory describing the first step of the 
process. 

One should find the solution at £ < and impose the boundary conditions (flUj) . Note, 
however, that the energy of our solution is fixed already. As for the initial oscillator excitation 
number z/, it does not change during the evolution in the region £ < 0. Thus, one may fix it 
at the matching time t = t\. One writes, 



"= \{y 2 + y 2 ) 



cos 2 j3 + sin 2 (3 sin 2 ^! cos (3) 



t=u 



This complex equation allows one to express t\ as 



sin^x cos/3) = —i — ^— — , (31) 
sin p 

where the choice of the sign is dictated by the condition in footnote [9j It is convenient to 
introduce notation t\ = iT\, T% 6 K. 

In order to find the suppression exponent fpiy), one needs to evaluate the parameters 
T(z/), 6{u). At £ < the solution has the form, 

x-(t)= po (t - iT/2) + x' , (32a) 
y _(t) = ae~ lt + a*e T+9+it , (32b) 

where the definitions (fl3|) . (fbil) have been taken into account explicitly, so that po, x' G R. 
One evaluates po, x' Q , a, T, 6 by matching the coordinates x±, y± and their first derivatives 
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Figure 3: The suppression exponent fp{y) for the waveguide ( 1201) ; (3 = tt/3. 



x±, y± at t = iT\\ this yields 



x' = 



Po 



nn T 

T x 

2 



= \/2(l - u) 
1 — vj cos 2 (3 



I V T+e 

a = i \ I - e 2 



sh 7\ 



\/ cos 2 (3 — v 
I — v \ 2 J sin (3 \fv 

The last two equations, together with Eq. fl25|) . define the function fpiy). 



2 f V^cr - f 

arcsh 



cos /3 



sin /3 



z/ cos /3 arcsh ,^ cr ^ - y/ (v CT - v) (1 - v) \ ; 



this finction is plotted in Fig. [3j One observes that at v — > z/ cr the quantities T%, T, 8, fp 
tend to zero, and the complex trajectory tends to the classically allowed critical solution, cf. 
Eqs. (EHD, 

Po — > \/2 sin (3 , a — > — p cos /3 . 

V2 



At z/ = one has, 



/fl(0) = -2 + arcth (cos (3) . 

cos p 



(33) 



To summarize, we obtained the suppression exponent for the reflection of a particle in the 
simplest waveguide with one sharp turn. 
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Figure 4: The equipotential contour w 2 (x, y) = 2N' for the waveguide fl35l) and the trajectory 
of the critical solution with energy N'/ cos 2 (3 > Eb- The matching points C, C are shown 
by the thick black dots. 

4 The model with two turns 
4.1 Introducing the system 

In the model of the previous section the suppression exponent was proportional to energy 
because of the coordinate rescaling symmetry (l23il . Now, we are going to demonstrate that 
small violation of this symmetry results in highly non-trivial graph for F(E). 

One introduces a second turn into the waveguide, see Fig. HI We want to consider this 
turn as a small perturbation, so, we assume its angle a to be smaller than (3. It is convenient 
to introduce two additional coordinate systems, x', y' and £, rj, bound to the central and 
rightmost parts of the waveguide respectively. They are related to the original coordinate 
system x, y as follows, 

x' \ / cos a sma \ ( x \ I £ \ / cos (3 — sin (3 

y' J \ — sin a cos a J \ y J \ r\ J lsin/5 cos f3 

' (34) 

Note that the origin of the coordinate system £, t] is shifted by the distance L. The waveguide 
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function is 

w = e(- x ')6(-^)y + 0{-£)6{x')y' cos a + 0(0*7 cos a cos (3 ; (35) 

it consists of three pieces glued together continuously at x' = and £ = (lines A'S' and 
AB in Fig. H] respectively). At t — > — oo the particle comes flying from the asymptotic region 
x' < 0, where w = y. In the intermediate region x' > 0, £ < the particle moves in the x' 
direction oscillating along the y' coordinate with the frequency cos a. Finally, in the region 
£ > its motion is free in the coordinates £, rj; the frequency of 77-oscillations is cos a cos f3. 

The model (135]) no longer possesses the symmetry (123]) : rescaling of coordinates changes 
the length L of the central part of the waveguide. In what follows it is convenient to work 
in terms of the rescaled dynamical variables, 

x = x/L, y = y/L. 

In new terms the parameter L disappears from the classical equations of motion, entering 
the theory through the overall coefficient L 2 in front of the action. The initial-state quantum 
numbers are also proportional to L 2 , 

E = L 2 E , N = L 2 N . (36) 

Thus, the conditions ( 1221) for the validity of the semiclassical approximation are satisfied in 
the limit 

L — » 00 , E, N — fixed . 
The suppression exponent takes the form 

F(E,N) = L 2 F(E,N) . (37) 

To simplify notations, we omit tildes over the rescaled quantities in the rest of this sec- 
tion. Rescaling back to the physical units can be easily performed in the final formulae by 
implementing Eqs. (136]) . ( |37j) . 

4.2 Classical evolution 

Let us begin this subsection by demonstrating that there are no stable classical solutions 
localized in the region between the turns. This is important for the determination of the 
tunneling probability, since such stable solutions could lead to exponential resonances in the 
tunneling amplitude. The argument proceeds as follows. Any trajectory which is localized in 
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the intermediate region should reflect from the line AB infinitely many times. Each reflection 
involves touching the unstable orbit living at the line AB. This implies that the trajectory 
itself is unstable. 

We proceed by determining the region of initial data E, N, which correspond to the 
classical reflections. [For brevity we will refer to this region as the "classically allowed 
region" , as opposed to the "classically forbidden region" where reflections occur only at the 
quantum mechanical level. We stress that these are the regions in the plane of quantum 
numbers E, N.] Let us search for the critical classical solutions which correspond to the 
smallest initial oscillator number iV = N cr (E) at given energy E. As in the previous section, 
one finds that the particle must get stuck at the lini^l AB for some time in order to reflect 
back. Let us first make an assumption inspired by the study of the one-turn model that the 
critical solutions touch the line AB at their maximum inclination point (point C in Fig. H]). 
We will see shortly that this is true only at energies above a certain value Eb, see Eq. ( 1501 . 
Still, the analysis based on the above assumption enables one to catch the qualitative features 
of the critical line N = N CT (E). Besides, the analysis is considerably simplified in this case; 
we postpone the accurate study until the end of this subsection. Keeping in mind the above 
remarks, one writes for the solution in the intermediate region, 



Before entering the intermediate region, the particle crosses the line A'B' (point C in Fig. [4]). 
The initial oscillator number N is most conveniently calculated at the moment 




(38a) 




(38b) 



cos a 



1 



t = t = 




of crossing. Using the relations (IMj) one obtains, 




(39) 



and thus 



N a (E) = E — -x 2 CT (to) = E — E sin (3 cos a — cos (3 sin a cos 



cos a 



E > E 



V2Esin(3 



B ■ 



(40) 



17 



0.2 



0.15 - 




^ ' ^ ' ' ' ' ' ' 1 

0.1 E B 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

E/oc 2 

Figure 5: The boundary N = N cr (E) of the classically allowed region at E > Eb for the 
waveguide model fl35l) ; (3 = ir/3, a = 7r/30. The region of the classically allowed initial data 
lies above this boundary. The empty circles correspond to the energies E = E n , where the 
curve N = N cr (E) touches its lower envelope N = Ecos 2 (/3 + a). 

As an example, we show in Fig. the region of the classically allowed initial data for (3 = 7r/3, 
a = 7r/30. One observes that the function N cr (E) oscillates between two linear envelopes, 
Ecos 2 (P + a) and E cos 2 \(3 — a); the period of oscillations decreases as E — > 0. Moreover, 
the curve N cr (E) has a number of minima at the points E = E™. This means that the 
energies E = E™ are optimal for reflection: in the vicinity of any point E = E™, N = 
N cr (E^) reflections become exponentially suppressed independently of whether the energy 
gets increased or decreased. This feature is particularly pronounced in the case ot + (3 = tt/2, 
when the lower envelope coincides with the line N = 0. Then, the classical reflections (i.e. 
reflections with the probability of order 1) at N — are possible only in the vicinities of the 
points 

8vr 2 (n- 1/2) 2 ' 

This is the case we used in Introduction to illustrate the effect. 

The minima E = E™ exist at other values of the parameters as well. For instance, let 

13 We do not consider reflections from the line A'B'. They disappear at larger values of N than reflections 
from the line AB if a is small enough. 

18 



us find the positions of these minima in the case a < 1. One differentiates Eq. (140!) with 
respect to energy and obtains, 

E " = E " { 1 - ^wr sin ( 2^^-1/2) ) + °< a2) } • (41 > 

where 

En = 87r2(n- 1/2)2 sin 2 /3 ( 42 ) 
are the points where the curve N = N cr (E) touches its lower envelope. The argument of 
arcsine in Eq. (|41|) should be smaller than one, so, the minima E" exist only at large enough 
n. 



n > n 



ctg/3 1 



+ 1 , (43) 



^ 2wa 2_ 
where [•] stands for the integer part. 

Let us make several comments. First, note that no ~ 0(1 /a), consequently, all the 
optimal points E™ lie in the region of small energies E ~ l/^o ~ 0(a 2 ). Second, as we 
pointed out before, the formula fj40|) for the function N CT (E) holds at E > Eg. Comparing 
the expressions flUD , (H3|) and (150|) . one observes that E no > E B if tg /? > 1. So, there 
does exist a range of energies where the non-monotonic behavior of the function N cr (E) can 
be inferred from the formula pUl) . In fact, the conclusion about the existence of the local 
minima of N cr (E), as well as the expressions ( HIT) . ( f42l) . ( 1431) determining their positions, 
remain valid also at E < Eb- This follows from the rigorous analysis of the boundary of the 
classically allowed region to which we turn now. The reader who is more interested in the 
tunneling processes may skip this part and proceed directly to subsection 14.31 

Now, we do not appeal to the Ansatz (!38j) . Instead, we start with the general solution 
in the intermediate region, 

x' =p' Q (t- t ) , (44a) 
y = A' Q sin [(£ — to) cos a + <p'] . (44b) 

It is convenient to parametrize it by the total energy E = p' 2 /2 + cos 2 aA 2 /2 and the 
"inclination" 7 defined by the relation 

Po/ A 'o = tg7cosa . 

Expressions ( [441) take the following form, 

x = V2E (t- t ) sin 7 , (45a) 
cos 'y 

y' = V2E sin \(t - t ) cos a + tp'\ . (45b) 

cos a 
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The constants to and ip' are fixed by demanding the trajectory (1431) to reflect classically from 
the second turn, i.e. touch the line £ = at t = 0, 



[x 



1) cos (3 — y' sin /3 
These conditions imply, 







t=o 



ctg/? 



t=0 




arccos 



tg7 
tg/3 



(46a) 



(46b) 



One sees that the classical reflections are possible only at 7 G [0; /3]; the boundary value 
7 = (3 reproduces the solution (138|) . 

In order to find N CT (E), one should minimize the value of the incoming oscillator exci- 
tation number with respect to 7 at fixed E. At t = to, when the particle crosses the first 
turn, 

Po = x(t ) = V / 2~E(cosa sin 7 — sin a cos 7 cos <p') . (47) 

Since N = E — one can maximize the value of the translatory momentum p instead 
of minimizing N(^y). Formula (T3TJ|) represents the value 7 = /? lying at the boundary of the 
accessible 7-domain; this value should be compared to Po(7) taken at local maxima. 

Let us consider the case a C 1. At large enough energies, E ~ 1, Eq. (|17|) is dominated 
by the first term, which grows with 7, so that the maximum of Po(7) is indeed achieved at 
7 = /3. At small energies, however, the second term in Eq. (14T1) becomes essential because of 
the quickly oscillating cosy/ multiplier: the frequency of cosy?' oscillations grows as E —>■ 0, 
and at E ~ a 2 , in spite of the small magnitude proportional to sin a, the second term 
produces the sequence of local maxima of the function poi'j)- 

One expects the parameters of the trajectory at small a not to be very different from the 
ones at a = (the latter case was considered in Sec. [3]). So, we write, 

T = (3- 5 1 , 

where < £7 <C 1. Expanding the expressions (|46l) . (j4"71) and taking into account that 
E ^ a 2 one obtains, 

1 



;i + 5 7 ctg/3) 



V2^sin^ 
2_E(sin /? — 57 cos (3 — a cos /3 cos <//) 



(48a) 
(48b) 
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Now, the local maxima of the initial translatory momentum can be obtained explicitly by 
differentiating Eqs. (j48p with respect to 8j. One finds the sequence of them, 



5 ln = -tgp + V2E 



sin 2 j3 
cos (5 



2irn 



71 



arcsm 



fZE sin 2 
a cos P 



(49) 



Only the maxima with £7^ > should be taken into account. The local maxima exist when 

a 2 cos 2 P 



E <J Eg = 



2 sin 4 P 



(50) 



Substituting Eq. (j4"9l into the expressions flSD , one evaluates the values of po at the local 
maxima, 



Po n (E) =2\f2E sin p - 2E sin 2 p 



2nn — Ti — arcsin 



fZE sin 2 P 
a cos P 



+ aV2EcosP\ 1 - 



2E sin 4 P 



a 2 cos 2 P 



The graphs N n (E) = E — p\ n (E)j2 are shown in Fig. [6] for the case P = tt/3, a = tt/30. 
Each graph is plotted for the energy range E > Ea u restricted by the condition 5^ n > 0. 
They are presented together with the curve given by the formula (j4"0|) . By definition, the 
critical solution corresponds to the lowest of these graphs. Clearly, for each "local" curve 
representing the n-th local minimum of ^(7) there is a range of energies E^ n < E < Eb„ 
where it lies lower than the "global" curve flUl) . This means that the parameter 7 of the 
critical solution changes discontinuously across the points E = Eb„- Correspondingly, the 
curve N cr (E) has a break at these points. On the other hand, the function N cr: (E) is smooth 
at the points A n as the "local" graphs end up exactly at £7 = 0, where the parameters of 
the n-th "local" solution coincide with the ones of the "global" solution. 

To summarize, we have observed that the boundary of the classically allowed region is 
given by a collection of many branches of classical solutions, each branch being relevant in 
its own energy interval. We will see that a similar branch structure is present in the complex 
trajectories describing over-barrier reflections in the classically forbidden region of E, N. 



4.3 Classically forbidden reflections 

In this subsection we demonstrate that the suppression exponent F(E, N) viewed as a 
function of energy at fixed TV exhibits oscillations deep inside the classically forbidden region 
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Figure 6: The graphs N n (E) corresponding to the local minima of the function iV(7) (dashed 
lines) plotted together with the "global" curve, Eq. (1401 (solid line); (3 = ir/3, a = tt/30. 
The critical curve N = N CT (E) is obtained by taking the minimum among all the graphs. 



of initial data. This result comes without surprise if one takes into account the non-monotonic 
behavior of the boundary N cr (E) of the classically allowed region. Indeed, the curve N = 
N cr (E) coincides with the line F(E,N) = 0. One has, 



dN cr 
dE 



d F F 



d N F 



N=N CT (E) 



so that 



dF 
dE 



We conclude that the points E = E^ are the local minima of the function F(E) at fixed 
N = N™. It is natural to expect that such local minima of F(E) exist at other values of 
N as well. To illustrate this fact explicitly, we study the complex trajectories, solutions to 

Eqs. ®, (HUD, (H2D. 

Following the tactics of the previous section, we find solutions in three separate regions: 
initial region x' < 0, final region £ > 0, and the intermediate region x' > 0, £ < 0. These 
solutions, together with their first derivatives, should be glued at t = to, when the complex 
trajectory crosses the line x' = 0, and at t — t±, when £ = 0. Besides, we are looking for the 
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tunneling solution which ends up oscillating along the line AB, see Fig. HJ As discussed in 
Sec. Elthis assumes existence of the second step of the process: classical decay of the unstable 
orbit living at £ = 0; the latter decay is described by a real trajectory 14 ! going to x — > — oo 
at t — > +00. 

The solution in the final region £ > is (cf. Eqs. ( 13T)j) '). 

£+(*) = , (51a) 
V2E 

ri+(t) — shift cos a cos /3) , (51b) 

cos a cos p 

where we used the time translation invariance (118p to fix the final oscillator phase (p v = 0. 
In the intermediate region x' > 0, £ < one writes, 

x'(i) = p' t + x' Q , (52a) 
y '(t) = a'e- itcosa + a'e ttcosa . (52b) 

Note that the final solution (j5ip does not contain free parameters; thus, the matching of x', 
x', y', y' at t = ti enables one to express all the parameters in Eqs. ( 1521) in terms of one 
complex variable ti, 

Pq = V2E sin /3 cos <j>\ , (53a) 

s'n = 1 + v 2_E [sin 0i — 0i cos 0i] , (53b) 

cos a 

a' = V^Z? e ^ l/cos/3 [sin 0i + i cos /? cos 0i] , (53c) 
cos a 

_, = V5ZH e _^ l/cos/3 r _ icos/3cos 0l ( 53d ) 

cos a 

where we introduced 0i = t\ cos a cos (3. 

As the energy of the solution has been fixed already, the only remaining initial condition 
involves initial oscillator excitation number at x' < 0, see Eqs. f 1 101) . It is convenient 
to impose this condition at the matching point t = to. One recalls the definition of the 
matching time to, 

p' t Q + x' = , 



14 One wonders why this trajectory does not reflect from the turn A' B' on its way back. This concern is 
removed by the observation that the trajectory produced in the decay of the unstable orbit is not unique: 
in appendix we show that the decay can occur at any point of the segment AC giving rise to a whole bunch 
of potential decay trajectories. Most of these trajectories pass through the turn A'B' without reflection. 
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which, after taking into account the expressions (I53al) . (I53bl) . leads to the following equation, 

cos a sin 0i . „ . 

+ -j - cos 0! A0 = , (54) 



y/2Esm(3 cos/? 
where A<p = cosa(ti — t Q ). At t = t one has, 

x(to) = p'o cos a — y(t ) sin a = y/2(E — N) , 

and thus 

= ctg a sin (3 cos <p\ — sin 0i sin A0 — cos (3 cos 0i cos A0 . (55) 

sin a 

As before, v = N/E. 

Two complex equations ( |54l) . ( !55l) determine the matching times to ; and, consequently, 
the complex trajectory. Although these equations cannot be solved explicitly, they can 
be simplified in the case a <C 1, which we consider from now on. For concreteness, we 
study reflections at iV = 0. It is important to keep in mind that in the region of interest 
E ~ E™ ~ 0(a 2 ); thus, one should regard all the momenta p and oscillator amplitudes a, 
a, as the quantities of order 0(a). At the same time, for the distances along the waveguide 
one has x ~ 0(1), so that the real parts of time intervals may be parametrically large, 
Ret ~ x/p ~ O(l/o). 

Further on, it will be convenient to work in terms of real variables, so, we represent <pi 
and A0 as 

(pi = cosacos/5(ri + iTi) , A(f> = cosa(r + iAT) . 

Note that r and AT are the real and imaginary parts of the time interval t\ — to which the 
particle spends in the intermediate region. Now, equation fl54l) enables one to express 



2E sin /Mi(Ti cos (5) 



+ 0(a), (56) 



T COS (3 



cos P 



ATcth(Ti cos (3) 



+ 0(a 3 ). (57) 



Note that T\ ~ 0(a), r ~ 0(1 /a). Then, the real part of Eq. (|55|) implies that 

ch(T lC os/5) = [l + actg/3cosre AT ] + 0(a 2 ) . (58) 
sm l) 

While deriving this formula we imposed T\ < which follows from the requirement that 
in the limit a — > equation fl3T|) should be recovered; besides, we assumed e AT ~ 0(1). 
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Substituting Eq. (!58|) into Eq. (156!) and the imaginary part of Eq. (I5"5l) . we obtain the final 
set of equations, 



1 - tV2E = a ctg/3 cos re AT + 0(a 2 ) , (59a) 
(1 + AT)e" AT = a ctgpr sin r + 0(a) . (59b) 

These two nonlinear equations, still, cannot be solved explicitly. Nevertheless, one can get 
a pretty accurate idea about the structure of their solutions. 

Before proceeding to the analysis of the above equations, let us derive a convenient 
expression for the suppression exponent Fq(E) = F(E,N = 0). Note that on general 
grounds one expects to obtain an expression of the form, 

F (E)=E(f p (0)+O(a)), 

where fpiO) is given by Eq. ( 1331) . We are interested in the 0(a) correction in this expression, 
so, one must be careful to keep track of the subleading terms during the derivation. 

Making use of the equations of motion, one obtains for the incomplete action (fl6|) of the 

system, 

2 Im S = Im p' = \JlE sin (3 Im(cos <pi) . 
Substitution of Eqs. (15"6"j) . (15T1) . (I5"5j) into this formula yields 

2 Im S = 2E I —1 — AT — a ctg j3 cos re AT ( 1 + — L- + 2AT ] + 0(a 2 
{ \ cos A p J 



For the parameter T one has (see Eqs. ([T 

T = - 2 -^ = - 2Im(x(to) - p ° to) = 2( Tl - AT) + .flsmalmy'ito) , (60) 
Po Po V E 

where in the last equality we used Eqs. (IM]) and x'(t Q ) = 0. The quantity hay' (to) is 
evaluated by using Eqs. (I52bj) . (|53|) and (1581) ; one finds, 

lmy'(t ) = -v / 2E(ctg (3 cos re AT + 0(a)) . 

Substituting everything into the formula (Tl5|) . we obtain, 

F Q (E) = E{jp(Q) -4actg/3cosrATe AT + 0(a 2 )) . (61) 

This expression implies that determination of the 0(a) correction to the suppression expo- 
nent involves finding r, AT with 0(l)-accuracy. This is precisely the level of accuracy of 
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Eqs. (I59p . Below we will also need the following formulae, which can be easily obtained by 
using T = -gf and Eq. (IBDI), 



dF 



dE 

_d_ (Fo 
dE \E 



f f3 (0) + 2(AT + l)+O(a), 
2(AT + l + 0(a)) 



E 



(62) 
(63) 



Note that, though the suppression exponent differs from that in the one-turn case only by 
0(a) correction, its derivative gets modified in the zeroth order in a. 

Now, we are ready to analyze Eqs. (|59|) . One begins by solving Eq. (159bl) graphically, 
see Fig. [7J The important property of this equation is as follows. One notices that the l.h.s. 
of Eq. f!59bj) is always smaller than 1, the maximum being achieved at AT = 0. Therefore, 
the solutions to this equation are confined to the bands 

tg/3 



r sin t < 



a 



This corresponds to 



r G [0; 27r(ni - 1) + 6r n 



or 



t G [2nn — n — <5t„; 2irn + 5r n ] , n > n\ (64) 



where 



5r„ 



arcsm 



tg(3 



tgff 1 
2na 2 



2na(n - 1/2) 
1 , 



0(a) 



(65) 



with [•] in the last formula standing for the integer part. The forbidden bands, where 
rsinr > tg/3/a, are marked in Fig. [7J by yellow shading. The property f)64p introduces a 
topological classification of the solutions r, AT to Eqs. fl59|) . Namely, these solutions fall 
into a set of continuous branches: the "local" branches T n (E), AT n (E) living inside the 
strips r G [27m — tt — 5r n ; 2irn + <5r n ], n > ni, and the "global" branch r g (E), AT g (E) 
inhabiting the very first band r G [0; 2n(ni — 1) + 5r ni ]. As follows from the definition of 
r, the topological number n counts the number of y'-oscillations during the evolution in the 
intermediate region. 

Let us consider the "global" branch. From Eqs. (159]) one has, 



2vr(ni - 1) + 0(a\na) 
. 



AT 3 
AT, 



ln(tg/5/a) 
-1 , 



E 
E 



0. 

+oo 
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Figure 7: Curves representing solutions to Eq. fl59b|) ; j3 = ir/3, a = 7r/30. 



By inspection of Fig. [7] one can work out the qualitative behavior of the functions r g (E), 
ATg(E). Alternatively, these functions can be found numerically. They are plotted in Fig. [8] 
for the case j3 = tt/3, a = tt/30 (the curves marked with "g"). One observes that at 
high enough energies the function AT g (E) exhibits oscillations around the line AT = — 1. 
According to the formula (163|) this means that the function F Q (E)/E is non-monotonic, it 
attains local minima at the points 



E' 



1 



Moreover, if 



n > n' = 



rr 2 (n- 1/2) 2 

"tg/3 



[1 + 2ae- 1 ctgf3 + Oi 



a 



Ana 



/^(O)exp 1 



2/2 



+ 1 



(66) 



(67) 



there exist E° n = E' n (l + 0(a)), such that AT(E°) = -1 - fp(0)/2. Then, according to 
Eq. (1621) the points E° are the "optimal" energies corresponding to the local minima of the 
suppression exponent F (E). 

At low energies the function AT g (E) ceases to oscillate and becomes large and positive. 
According to Eq. ( 1621) this means that the suppression exponent Fo^(E) of the "global" 
solution becomes negative at low energies 3, see Fig. [9j This is a clear signal that the 
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It is worth mentioning that Eqs. (|59[) and the expression (|61[) for the suppression exponent become 
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Figure 9: The suppression exponent F Q (E) for the "global" and first "local" (n = 4) branches; 
(3 = 7r/3, a = 7r/30. The vicinity of intersection of the graphs is enlarged in the upper right 
corner. 



"global" solution becomes unphysical at these energies and its contribution to the reflection 
probability should be discarded: negative suppression exponent contradicts the unitarity 
requirement! 16 !. V < 1. One is forced to conclude that at low energies reflection is described 
by the "local" solutions. Let us study them in detail. 
For the n-th branch one obtains, 

r n -> 2vrn + O(alna) , AT n -> In(tg/3/a) , E -> , 

r n — > 2™ — 7r , AT 9 — > +oo , — > +oo . 

From Fig. [7] one learns that the n-th solution passes through the points 

AT n = — 1 , t = 27tn or t — 2nn — n . (68) 



inapplicable at large AT: the assumption e AT ~ 0(1) which was used in the derivation of these equations 
gets violated. Nevertheless, by analyzing the full equations (|54|) . (j55]l one can show that dF ^ g /dE = —T g 
is large and positive at E — * 0. This is sufficient for concluding that Fo >g (E) is negative in the low-energy 
domain. 

16 Another indication that the "global" solution is unphysical at small E is that the function r g (E) is 
bounded from above. Indeed, r is the time interval the particle spends in the intermediate part of the 
waveguide, one expects it to tend to infinity as E — > for a physically relevant solution. 
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Thus, each curve AT n (E) has one sharp dip, its minimum is smaller than —1, see Fig. [HJ As 
in the case with the "global" branch, the points (1681) represent the extrema of the functions 
Fq jU (E) / E; the positions of the local minima are again given by Eq. (l66j) . 

Making use of Eq. (IBTj) . we find that the suppressions F 0yn (E) of the "local" branches are 
large and positive at high energies. Hence, these solutions give subdominant contributions 
to the reflection probability at such E as compared to the "global" solution. As energy 
decreases, Fo >n (E) also decreases, then makes one oscillation and drops to negative values at 
small E. The latter property means that each "local" branch becomes unphysical at small 
enough energies. The suppression exponent of the first "local" branch (corresponding to 
n = 4 in the case (3 — 7r/3, a — n/30) is presented in Fig. ED 

An alert reader may have already guessed that we have met here the typical Stokes 
phenomenon [21]. In fact, the Stokes phenomenon is specific to the situations where some 
integral (e.g., the path integral <^ in our case) is evaluated by the saddle-point method. 
Essentially, it means the following: as one gradually changes the parameters of the integral 
in question, a given saddle point may become non-contributing after the values of these 
parameters cross a certain curve drawn in the parameter space, the Stokes line. Since the 
result of the computation should be continuous, this phenomenon occurs only for subdomi- 
nant saddle points (saddle-point trajectories in our case). Unfortunately, apart from several 
heuristic conjectures (2T], [12] , sometimes rather suggestive [13], there is presently no general 
method of dealing with the Stokes phenomenon in the semiclassical calculations. However, 
in the situation encountered above it suffices to use the simplest logic lying at the heart of 



all other approache 

When gathering the final result for the suppression exponent, we follow two guidelines. 
First, it is clear that, as energy decreases, each branch becomes unphysical before F 0;n (E) 
crosses zero. On the other hand, at high energies one should pick up the branch corresponding 
to the smallest value of the suppression exponent. Looking at Fig. [9j one notes that the 
curves Fq^{E), Fq^{E) have two intersections, A and B. At E > Eb one chooses the 
"global" branch. In the region Ea < E < Eb we switch to the first "local" branch, because 
in this region Fqa(E) < Fo j9 (E). Naively, at E = Ea one should jump back to the "global" 
branch; however, in order to preserve unitarity at small energies, we suppose that somewhere 
in between the points B and A the "global" branch becomes non-contributing, so that one 
should stay at the "local" branch at E < Ea- Similarly, the adjacent "local" branches have 

17 The simplification in the present case is related to the fact that we concentrate on the dominant semi- 
classical contribution, leaving aside the subdominant ones. 
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two intersections; as the energy decreases, we switch from n-th branch to n + 1-th at the 
first intersection, and stay there until the intersection with the n + 2-th branch. Overall, one 
obtains the graph for the suppression exponent plotted in Fig. [TUJ The suppression exponent 

0.08 i 1 1 1 1 1 




1 '- ' ' - i 2 i i i 

0.05 0.1 0.15 0.2 0.25 



E/cc 2 

Figure 10: The final result for the suppression exponent F (E) in the region of small energies; 
(3 = 7r/3, a = 7r/30. The points where different branches merge are shown with thick black 
dots. 

oscillates between two linear envelopes, F = E(fp{Qi) ± 4e~ 1 a ctg/3); oscillations pile up in 
the region of low energies. The reflection process is optimal in the vicinities of the minima 
of the function F (E). 

5 Discussion 

By considering a class of two-dimensional waveguide models, we have demonstrated explicitly 
that the probability of over-barrier reflection can be non-monotonic function of energy. The 
origin of the effect lies in the classical dynamics: the parameters of the complex trajectory 
describing over-barrier reflection change quasi-periodically as the energy gets decreased. 
This results in the oscillatory behavior of the suppression exponent. Reflection occurs with 
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exponentially larger probability in the vicinities of "optimal" energies (local minima of the 
suppression exponent) while being highly suppressed in between. 

Our results are obtained for a fairly specific class of waveguides, namely, the ones with 
very sharp turns. However, the qualitative features observed in this paper should be valid 
for quite general waveguide models: a classical particle with high energy feels any large-scale 
turn of the waveguide as a sharp one^l; if two turns are separated by a long interval of free 
motion, one arrives to the model ( |35|) . We remark that the phenomenon of optimal tunneling 
has been observed also in numerical investigation of a smooth waveguide, see Ref. |14j . 

The branch structure of solutions observed in the region of small energies is interesting 
from the mathematical point of view. We have shown that there exists an infinite sequence of 
complex trajectories marked by the topological number n. Each branch produces physically 
consistent result for the suppression exponent in some energy interval; outside of this interval 
the n-th branch would correspond either to highly suppressed transitions (high energies) or 
to violation of unitarity (low energies). We collected the final graph for the suppression 
exponent basing on the empirical considerations, which hardly may be acknowledged as 
satisfactory. Our study clearly shows that the method of complex trajectories should be 
equipped with a convenient rule to pick up the physical trajectory among the discrete set of 
solutions to the boundary value problem (jHJ), (ITU]) . (1T2"]) (in other words, the method to deal 
with the Stokes phenomenon). Presently, such a rule is absent. 

We note that the described physical phenomenon of optimal tunneling is present inde- 
pendently of the way the branches of solutions are glued together. The result at relatively 
high energies is given by the "global" branch, which displays a large number of local minima 
if n' > m, see Eqs. ( )67|) . (1651) . This is the case for the illustrative example considered 
throughout this paper, see Fig. [91 

As a final remark, we point out some open issues. We have calculated the suppression 
exponent of reflection using the sharp-turn approximation. It would be instructive to extend 
our analysis by finding corrections due to the finite turn widths. The motivation is twofold. 
First, the analysis performed in appendix implies existence of a rich variety of distinctive 
semiclassical solutions contributing almost equally into the reflection probability. This fea- 
ture might be a manifestation of chaos [7] which is present in our system but hidden by the 
sharp-turn approximation. [Note that chaos is inherent in a very similar waveguide model 



18 More precisely, one should compare the width b of the turn to the quantity 77^, where po is the 
translatory momentum of the particle and u> stands for the frequency of transverse oscillations; if b <C wpo , 
one is in the class of models with sharp turns. 



32 



with smooth potential, see Ref . [H] .] Clearly, the structure of solutions in the vicinities of 
the turns is worth further investigation. 

Second, it was proposed recently in Refs. [HI [11] that the process of dynamical tunneling 
in quantum systems with multiple degrees of freedom (including field theoretical models, 
see Refs. [TJJ]) can proceed differently from the ordinary case of one-dimensional tunnel- 
ing. Namely, classically unstable state can be created during the process; this state decays 
subsequently into the final asymptotic region. The analysis performed in the present paper 
naturally conforms with this tunneling mechanism: all our complex trajectories are matched 
with the unstable orbit living at the turn. Still, the sharp-turn approximation does not allow 
to distinguish between the truly unstable trajectories staying at the turn forever and those 
which reflect from the turn in a finite time. To decide whether the tunneling mechanism of 
Refs. [91 [11] is indeed realized in our model one needs to go beyond the sharp-turn approx- 
imation. Then, the candidate for the "mediator" unstable state is the "excited sphaleron" , 
the solution considered in the appendix. Presumably, in our model one can answer analyti- 
cally to the question of whether or not the "excited sphaleron" acts as an intermediate state 
of the tunneling process. This study is quite beyond the scope of the present paper and we 
leave it for future investigations. 
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A Classical motion near the turn 

In this appendix we analyze the motion of the particle near the sharp turn of the waveguide 
( |20|) at nonzero smoothening of the turn, see, e.g., Eq. ( 12TI) . We suppose that in the small 
vicinity of the turn the function u>(£, if) can be represented in the form 

w(£,Tj) = caaP(r l -bo{£/b)) , (69) 
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where v (ip) does not depend explicitly on b. Moreover, we consider the case when v(ijj) has 



a maximu 



v'(ipo) = . (70) 



Due to the property (I70j) one immediately obtains the exact periodic solution to the 
equations of motion ([§]), which we call "excited sphaleron" [S], 

£ sp = bip , i] sp = A v sin(t cos/3 + (p v ) + bv(ip ) . (71) 

We are going to show that this solution is unstable: a small perturbation above it grows 
with time and the particle flies away to either end of the waveguide. In particular, there are 
solutions that describe the decay of the sphaleron to £ — > — oo both at t — > ±oo. Clearly, 
such solutions correspond to reflections from the turn. 

In the vicinity of the sphaleron the trajectory of the particle can be represented in the 
form, 

S = W(t), r ] = r ]sp (t) + bp(t) , (72) 

where ip, p ~ 0(1). Writing down the classical equations of motion ([9]) in the leading order 
in b, one obtains, 

d 2 ip 4 



A v sin(2s)i/(V>) , (73) 



2 +4p = 4K^)-^ )], (74) 



ds 2 b 
ds 



where s — (t cos (3 + ^>^)j1. It is worth noting that the right hand side of Eqs. (1731) . (1731) are 
of different order in b. We will see that due to this difference p = in the leading order in b. 
Let us first consider the linear perturbations above the excited sphaleron, 

ip = ifj + 5ip , 5i/j < 1 . 

Equation (173j) can be linearized with respect to Sip leading to the Mathieu equation 

d 2 

-—5ip + 2qsm(2s)6ip = , 
ds^ 

with canonical parameter q = —Iv^A^jb > 0. As q ~ 0(1/6) ^> 1, one can apply the WKB 
formula, 

_ AcosW 

y/dW/ds ' 



19 For the smoothening the properties ((Ml), dZOJ) hold with v(ip) = ^||, Vo » 1-28. 
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where \A\ -C 1, and 



W = y/2q [ ds v / sin(2s / ) 



t/J 

Note that we have chosen the solution symmetric with respect to time reflections, 

^(tt/2 - a) = 8ip(s) . (76) 

At s G [0; 7r/2] the exponent W is real and the particle gets stuck at ip ~ ^o, oscillating 
around this point with high frequency dW/ds ~ 0(b~ l l 2 ). At s < the solution fl75l) grows 
exponentially, meaning that the particle flies away from the excited sphaleron, 

Acos( ^(0)-7 r/4) p|w(a) _ w(0) , 
v/|dW/ds| 

In what follows, we choose Acos(W(0) — 7r/4) < 0, so that 5^ < at s < 0. Let us denote 
by s\ < the point where dip becomes formally equal to —1, 

Acoa(W(0)-ir/4) mai) _ wm 



^\dW/ds 



In what follows we suppose that s\ ~ 0(1), hence, A is exponentially small. Then, in the 
vicinity of this point, \s — s\\ <C 1, one has, 

<fy; = -exp (V-2gsin(2si)(si - s)) = - exp ^4i;g^sin(2si) ^~y=p ) • ( 77 ) 

We notice that 5ip evolves from exponentially small values to 5ip ~ 0(1) during the charac- 
teristic time \s — s\\ ~ 

When Sip ~ 0(1) the linear approximation breaks down and one has to solve the nonlinear 
equation (1731 . Using s = s x + 0(\/b) one writes 

^ = ^sin(2 Sl KW . (78) 

This equation permits to draw a useful analogy with one-dimensional particle moving in 
the effective potential V e ff(i/)) = — 4b~ 1 A v sin(2si)v(ip) (see Fig. [TTj). This auxiliary particle 
starts in the region near the maximum of the potential at (s — Si)/yb — ► +oo with energy 
E ~ Knaa; and rolls down toward ip — > — oo at (s— si)/Vb — > — oo. In this limit t> (■?/>) — > ^ tg/3 
and the solution takes the form 

V> = Ci + C 2 (s - si) + 26 _1 A, sin(2si) tg/3 (s - si) 2 . 
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Figure 11: The effective potential for Eq. (178]) . 



Note that the coefficients C±, G% here are not independent: they are determined by the 
parameter s\ through matching of the solution with Eq. (177)) at (s — sx)j\fb — > +00. We do 
not need their explicit form, however. 

Let us argue that the function p remains small during the whole evolution of the particle 
in the vicinity of the sphaleron. Indeed, in the linear regime one has Sip <C 1 and the r.h.s. 
of Eq. (1731) is small. So, p does not get excited. On the other hand, the nonlinear evolution 
of ip proceeds in a short time interval As = 0(y/b); so, again, p is suppressed by some power 
of b. 

The trajectory ( 1721) found in the vicinity of the sphaleron should be matched at 

1 > \s - 8 l \ > Vb 

with the free solution in the asymptotic region £ < 0, see Eqs. (1261) . It is straightforward to 
check that matching can be performed up to the second order in (t — ti), which is consistent 
with our approximations. In this way one determines the free asymptotic solution which, up 
to corrections of order 0(b), coincides with the sinusoid coming from £ cxoati^— c» 
and touching the line £ = at t = t\. 

Now we recall that, by construction, the obtained solution is symmetric with respect to 
time reflections, 

£(s) = £(tt/2 - s) , V (s) = V (n/2 - s) . 

This means that it satisfies £ — * — 00 at t — * ±00. This solution describes reflection of the 
particle from the turn. 
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The reasoning presented in this appendix puts considerations of the main body of this 
paper on the firm ground: we have found the "smoothened" solutions which reflect classically 
from the turn, and in the limit b — > coincide with the free solutions of Sec. [3] touching the 
line £ = 0. 

It is worth mentioning that, apart from the reflected solution we have found, in the 
vicinity of any trajectory touching the line £ = there exists a rich variety of qualitatively 
different motions. First of all, one may successfully search for solutions which are odd with 
respect to time reflections (Eq. (|76l) with minus sign). Such solutions, though close to the 
reflected ones at t < 0, describe transmissions of the particle through the sharp turn into the 
asymptotic region £ — > +oo. Relaxing the time reflection symmetry, one can find solutions 
leaving the vicinity of the turn at any point 77 < 0, which is different, in general, from the 
starting point rj = rj(si). Yet another types of solutions are obtained in the case when the 
amplitude A of (^-oscillations at s G [0; 7r/2] is so small that Sip does not reach the values 
of order one during the time period s G [—tt/2; 0]. If the particle is still in the vicinity of the 
point ipQ at s = —tt/2, it remains for sure in this vicinity at s G [— 7r; — 7r/2], because the r.h.s. 
of Eq. (1731) is positive again. In this way one obtains solutions, which spend two, three, etc. 
sphaleron periods at ip ~ ipo before escaping into the asymptotic regions ip — *> ±00. In the 
leading order in b all these solutions correspond to the identical initial state, and (in the case 
of classically forbidden transitions) to the same value of the suppression exponent. However, 
an accurate study of the dynamics in the vicinity of the the sphaleron is generically required 
to obtain the correct value of the suppression exponent in the case b ~ 1, cf. Ref. p3]. 
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